Possible stimulation of nuclear a decay by superfluid helium 
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It is suggested that superfluid helium (condensate of 4 He atoms) may stimulate nuclear a decay 
in a situation when an a emitter moves through superfluid helium with fine-tuned velocity, so that 
the backward-emitted a particle is at rest in the laboratory frame. It is shown that the probability 
of stimulated a decay in this case may be sizable enough to be detected. 
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INTRODUCTION 

Possibilities for influencing nuclear reactions by mod- 
ification of the environment have been discussed many 
times. In particular, the highly developed field of muon 
catalysis belongs in this domain. The purpose of this let- 
ter is to draw attention to the possibility, in principle, of 
stimulating nuclear a decay by superfluid helium in some 
specific conditions. 

In a Bose system the amplitude of a particle transition 
to a state is proportional to ^/N + 1, where N is the state 
occupation number. The transition probability, hence, is 
proportional to N + 1 and consists of two contributions 
caused by stimulated and spontaneous transitions. The 
stimulated transitions may result in boson accumulation 
in a single quantum state as, e.g., it occurs for photons 
in a laser. Bosc-Einstein condensation (it was discovered 
[3, 0] in dilute gases of alkali bosonic atoms in 1995) has 
the same basic explanation. 

Indeed, let us show that the temperature of Bose- 
Einstein condensation, Tb = (27r/i 2 /m)(n/2.612) 2 / 3 (see, 
e -g-i |3() results from a competition between spontaneous 
and stimulated transitions. Let T be the temperature of 
a Bose gas, V be the volume, n be the gas density, and 
the population of the ground state (Bose condensate) be 
iV c ~ nV. Spontaneous transitions cause the scatter- 
ing bosons to populate the states belonging the phase 
space volume AnVp^/S (the momentum px is defined, 
e.g., by p\j2m — 3T/2, where m is the boson mass). 
On the other hand, stimulated (by the condensate) tran- 
sitions leave bosons in the ground state. Transitions to 
the ground state dominate, when: 



N c > 



47rVp| 
3(27rft) 3 



T < 



h 2 /27T 2 n^ 2/3 



V3 



To- (1) 



Thus, the Bose condensate survives at low temperatures 
T < To- It can be seen that To and Tb differ only by a 
factor of the order of unity. 

Notice that the Bose-Einstein condensation occurs 
when the thermal de Broglie wavelength Xt = 2ttH/pt 
surpasses the interatomic distance, because Eq. ([1} may 
be rewritten as: 



i\?r > 4tt/3. 



(2) 



It is recognized nowadays that the superfluidity of liq- 
uid helium (the system of bosonic atoms 4 He) is related to 
the Bose-Einstein condensation. It is seen, in particular, 
from the closeness of the temperature Tb = 3.1 K (calcu- 
lated for the density of liquid helium n = 2.2 • 10 22 cm -3 ) 
to the temperature T\ = 2.17 K of helium transition to 
superfluid state. However, owing to the interaction, the 
atoms of superfluid helium populate not only the low- 
est quantum state (as particles of an ideal Bose gas), 
but some set of low-lying quantum states. Both theoreti- 
cal calculations and neutron scattering experiments show 
that at T < 1 K approximately 10% of the 4 He atoms are 
in the lowest state (see, e.g., 0, H, El, 0] and references 
therein). Thus, the density of the atomic condensate can 
be taken to be at least n c = O.ln. 

Let us consider now some nuclear process implying 
emission of an a particle - nucleus of a 4 He atom (or 
4 He atom if the emitted a particle is accompanied by 
two electrons) with the energy varying from zero up to 
the value of the order of T\. For the process in super- 
fluid helium we would expect an enhancement of the 
emission rate due to stimulated emission. However, it 
is hard to specify such a process. Indeed, typical energy 
release in nuclear processes is of the order of 10 6 -10 7 eV, 
whereas T\ ~ 1CP 4 eV. Let X a — 2irh/p a be the de 
Broglie wavelength of an a particle with the momentum 
p a = \/2mE a and energy E a . For E a ~ 9 MeV we ob- 
tain A Q ~ 0.5 • 10~ 12 cm. Then the parameter ricA 3 in 
the left-hand side of Eq. © is very small: - 0.2 • 10~ 15 . 

Formally, it is possible to consider an a decay with the 
energy release ~ 10~ 4 eV. But according to the Geiger- 
Nuttall law, the lifetime of such an a emitter would be of 
the order of ~ 10 10 y. Obviously, no reasonable increase 
of the decay rate would lower this lifetime to an observ- 
able value. Nevertheless the situation is not hopeless due 
to two circumstances. 

1. After nuclear decay, an a particle moves with the 
energy E a ~ 9 MeV (the corresponding velocity is 
v a = ^2E a /m ~ 2 • 10 9 cm/s). Let us assume 
that a container with superfluid helium moves as 
a whole with the same velocity v a in the same di- 
rection. Then for the decaying nucleus in the vol- 
ume of helium, the a particle would be in the same 
quantum state with the momentum p a — mv a as 
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all the N c = n c V atoms (and their nuclei - a par- 
ticles) of the condensate. Therefore the probability 
of such a decay would be enhanced by a factor of 
N c . In practice, of course, it is easier to accelerate 
the decaying nucleus up to the velocity v a with re- 
spect to the superfluid helium. Then the a particle 
emission strictly backwards would be stimulated. 

2. The probability of spontaneous emission is propor- 
tional to the number of available states in the phase 
space. In the right-hand side of the first formula 
of Eqs. ([1]) we have taken into account all single- 
particle states with energy up to p\/2m, where pr 
is the maximum momentum. Just for this case the 
momentum pt (as well as the appropriate energy 
T) has turned out very small. But the decaying 
nucleus emits a particles into a very narrow region 
of the momentum space. Indeed, an uncertainty of 
the a particle energy is of the order of h/r (typ- 
ically much less than T\ ~ 1CP 4 eV), where r is 
the nuclear lifetime (e.g. h/r ~ 10~ 9 eV when 
r ~ 1CT 6 s). Thus, stimulated emission may come 
to light. 

Of course, there is an important question here. Is the 
condensate of 4 He atoms simultaneously the condensate 
of their point-like nuclei - a particles? To my mind, the 
answer is "yes" and it follows specifically from the re- 
sults of experimental searches of the condensate in liquid 
helium by means of neutron inelastic scattering (see de- 
tails, e.g., in Q). Indeed, neutrons interact just with the 
nuclei of 4 He atoms. On the other hand, the probability 
of two electrons being picked up by an a particle (such a 
possibility was already noted above) could be sizable be- 
cause v a is of the same scale as the velocities of electrons 
of the a emitter. Therefore estimates for stimulated a 
decay seem worthwhile even if the condensate is formed 
of 4 He atoms. 

Thus, below we consider the situation of an a emitter 
moving through a superfluid helium with the velocity v a , 
so that the backward-emitted a particle is at rest with 
respect to the condensate. We give an estimate for the 
probability of stimulated a decay and discuss the possi- 
bility of an experimental search for this phenomenon. 

PROBABILITY OF STIMULATED a DECAY 

Let us show that the probability of a decay really 
is determined by a number of available final states and 
this probability might be increased by a condensate. It 
is convenient to consider the decay caused by a time- 
independent perturbation V. Then the probability of 
transition from an initial state ipi with the energy Ei to 
a final state tp n with the energy E n during a certain time 
t is defined by (see, e.g., [H) 

P(i^n,t) = 4\(4, n \V\^)\ 2 sm 2 (A n t/2h)/A 2 n , (3) 



where A„ = E n — Ei. For a fixed t, the right-hand side 
of Eq. ([3]) represents a peak with the center at A„ = 
and the width AE ~ h/t. Thus, if t is not small, the 
final states belong to a narrow layer in the phase space. 

The total transition probability is determined by the 
sum over states within the limits of the peak or, other- 
wise, by the product of the right-hand side of ([3]) and the 
number AG(E) = p(E)AE of the final states (p(E) is the 
density of states). In practice one replaces the sum by 
the integral over E. Its calculation gives: Pif(t) — Wift, 
where the transition probability per unit time is defined 
by Fermi's rule: 



2tt 



\^f\V\^)\ 2 p{E } ). 



(4) 



Here Ef = Ei (because V does not depend on t). 

Now we notice that the initial wave function of the 
a particle, ipi, is localized inside the nucleus (r < 
r±). The final wave function, t/jf, is distributed in a 
macroscopic volume V and, therefore, inside the nu- 
cleus takes the form: tpf(r) = V'/MV 'D /V , where 
D = exp(-2 y / 2m(U(r) - E a ) dr/h) is the penetra- 
bility factor (ri and T2 are the turning points), while 
ip'Ar) is a function of the order of unity. Thus we obtain 
from Eq. ([4]) the usual expression for the probability of 
spontaneous a decay per unit time: 



TO 



; Anp 2 a dp a 
{2TthfdE a 



(5) 



The quantity 1/to implicitly includes the factor AG (the 
number of available final states), and To is of the order of 
the oscillation period of an a particle in the nucleus. 

Let us now estimate the probability of stimulated a 
decay. Let ^n(1,2 . . .N) be the wave function of the 
condensate of N 4 He atoms or their nuclei - a particles 
(a qualitative description of 4" at has been given by Feyn- 
man 0). The initial wave function of N + 1 particles, 
normalized to unity and symmetric with respect to per- 
mutation of any two particles, looks like 



*i(l,2...JV+l) = 



1 



[*jv(l,2...iV)V>i(iV+l) 



+ V N (l...N-l,N+l)ilj i (N) + ... 
+ * W (2...JV,JV+1)^(1)]. 

(6) 

Evidently, the matrix element over the final condensate 
state = ^n+i and ^>i given by Eq. ([6]) is enhanced 
by the factor of y/N + 1 in comparison with the single- 
particle matrix element 

(*/l£v"(*)l*<) = n/FTT (^ N+ i\V\^ N ^i). (7) 

i 

The probability of transition per unit time, hence, in- 
creases by the factor of TV + 1 . 
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It should be noted that a decay is a tunneling pro- 
cess. Thus, it is instructive to show independently that 
the barrier penetrability for bosons increases due to sym- 
metrization of the wave function. Before turning to this 
task, we note that this fact might be demonstrated exper- 
imentally by modern techniques, that allow the prepara- 
tion of Bose condensates as well as single atoms in double 
well traps (see, e.g., [13, EI III El)- To do this, an initial 



state of a single boson in one well and Bose condensate 
in the other well should be formed. 

Let us consider a one-dimensional double well poten- 
tial, U(x) — U{—x). The left and right wells are sepa- 
rated by a barrier. Thus, the single-particle wave func- 
tion, ip s (x), of the lowest state with the energy E s is 
symmetric, ip s (—x) = ip s (x), whereas the single-particle 
wave function, i\) a (x) , of the next state with the energy 
E a = E a + SE is antisymmetric, i.e. ip a {— x) = —%l) a {x). 
Let i/j s (x) and i/j a (x) be positive at x < 0. Then the 
time-dependent single-particle wave function, 

i>(x,t) = (^ s (x)e-^ h + Mx)e- tEat/h ) A/2, (8) 

describes the particle localized in the left well at t = 0. 
But at t = irh/SE the same function describes the parti- 
cle localized in the right well. In the other words, during 
the time irh/SE the particle passes through the barrier. 

Let N bosons now be in the right well at t = 0, whereas 
one more boson is localized in the left well. The time- 
dependent wave function of the whole system may be 
written as 

*(a;i . ..x N+ i,t) = 

(y s (xi---XN + i)e^ N+1 ^ t / h (9) 
+ (-1)" * (n . . . x N+1 ) e -W+D^M) /y /2 , 

where ^ s and are the products of symmetric and anti- 
symmetric single-particle functions, respectively. Assum- 
ing that at t = the coordinates xi, £2. . . xn are positive 
(i.e., N particles are in the right well), the function (J9j) 
is nonzero only if x^+i is negative (i.e. (N + l)th par- 
ticle is in the left well). However, at t = irh/(N + 1)5E 
with the same assumptions about x\, X2- ■ ■ Xn, the wave 
function © describes the (N + l)th particle in the right 
well. Thus, indeed, the time of transition through the 
barrier decreases by the factor of (iV + 1) (to describe ac- 
curately the initial state when N bosons are in the right 
well, while one more boson is in the left well, one needs 
to introduce additional terms into the right-hand side of 
Eq. ([9|), but the quoted terms are leading). 

Let us now estimate the ratio of the probabilities per 
unit time for stimulated, w s t, and spontaneous, w sp , 
a emission. The simplest assumption is the following: 
among AG = p(E a )T a of available final states, where the 
line width T a — h/r a is determined by the nuclear life- 
time r a , there is only one state for the backward-emitted 



Table I: The results of calculation of parameter b and prob- 
abilities P„ p and P s t of spontaneous and stimulated a decay, 
during the flying of the nucleus through the superfluid helium, 
for several polonium and astatine isotopes; the calculations 
are performed for n c = 0.22 ■ 





E a (MeV) 


T sp (s) 




b 




p 
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211 Po 


7.44 


0.74 




8.6 ■ 10 5 


1.4 


■ 10" 
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0.9 


• 10" 
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212p 


8.79 


4.3-10" 


7 


0.46 


2.3 


■ 10" 


2 


4.3 


■ 10" 


-3 


213p 


8.38 


0.6- 10" 


5 


6.6 


1.6 


■ 1(T 


3 


2.3 


• 10" 


-3 


214 Po 


7.69 


2.4 ■ 1(T 


4 


2.7- 10 2 


4.4 


■ 10" 


5 


4.9 


• 10" 


-4 


212 At 


7.64 


0.45 




5.2 • 10 s 2.3 


• 10" 
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1.2 


• 10" 
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213 At 


9.08 


1.8 ■ 10" 


7 


0.19 


5.3 


■ 10- 
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• 10- 
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214 At 


8.81 


0.8 • 10" 
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0.85 


1.2 


• 10- 
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3.9 


• 10- 
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215 At 


8.03 


1.4 ■ 10" 
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1.6 ■ 10 2 


0.7 


■ 10" 


4 


0.6 


■ 10" 
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a particle populated by A^ c condensate particles. Thus, 
we get: 



w st /w sp = A^/ AG. 



(10) 



Taking r Q = l/(w sp + w s t), we obtain a quadratic equa- 
tion for w s t- Its solution, 



Wst 



w 



(Vr+26- l)/2, 



depends on the dimensionless parameter, 
b = n c \ 2 a v a T sp , 



(11) 



(12) 



that is defined by the condensate density n c , the wave- 
length X a = 2irh/p a , the velocity v a — p a /m, and the 
nuclear lifetime r sp = 1/ w sp with respect to spontaneous 
a decay. 

When b ~ 1, the probability of stimulated a decay is 
comparable to the probability of spontaneous decay. But, 
curiously, it can be easily achieved: the small parameter 
n c \ a ~ 1CP 15 can be neutralized by the factor v a T sp /X a , 
that surpasses 10 15 for r sp > 10~ 6 s. Moreover, for rel- 
atively long-lived nuclei with lifetime r sp ^> 10~ 6 s, we 
get b 3> 1, and therefore w s t 3> w sp - However, the stim- 
ulated a decay may occur only during the short time 
t x = x/v a when the decaying nucleus moves through a 
layer of superfluid helium of thickness x. Thus in any 
case there is no reason to expect large absolute value for 
the probability P s t — w s tt x of stimulated decay (we intro- 
duce also P sp = w sp t x as the probability of spontaneous 
decay during the same time t x ). 



NUMERICAL ESTIMATES AND DISCUSSION 



The value b 
10~ 6 s. In this case P s t 
x ~ 10 cm. 

If b <S 1, then 

w st / w sp ^b/2 



1 corresponds to the lifetime r s 



P, 



sp 



10 for realistic value 



w st — n c X 2 a v a /2. (13) 
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Notice that w s t does not depend on the lifetime r sp (but 
herewith w s t -C w sp = l/r sp ). If, e.g., r sp ~ 10~ 7 s, then 
P sp ~ 10" 1 and P st - 1CT 2 for x ~ 10 cm. 

At last, the case T sp » 10 -6 s (6 > 1) corresponds to 



w st /w sp ~ ^st - y n c X%v a /2r sp . (14) 

If, e.g., r sp ~ 10~ 5 s, then P st ~ 10~ 2 for x ~ 10 cm 
(whereas P sp ~ 10~ 3 ). 

As an illustration, the results of calculations performed 
for four polonium isotopes (the daughter products are 
lead isotopes) and four astatine isotopes (the daughter 
products are bismuth isotopes) are presented in the Ta- 
ble U (the decay parameters are from [HI])- For each iso- 
tope the dimensionless parameter b (TT2]) has been calcu- 
lated for the known energy E a (and, hence, v a and A Q ) 
and lifetime r sp = T 1 / 2 /ln2 at n c = 0.22 • 10 22 cm~ 3 . 
The probabilities P sp — x/(T sp v a ) and P st = w st x/v a of 
spontaneous and stimulated a decay are obtained for the 
layer thickness x = 20 cm (w st is given by Eq. HP) ). 
We see that the probability P st belongs to the interval 
from 10~ 3 to 10~ 2 when the parameter b is comparable 
to unity (in agreement with rough estimates presented 
above). 

Up to now we discussed the situation with the nucleus 
moving through the superfluid helium with the fine-tuned 
velocity v a , so that the backward-emitted a particle is at 
rest in the laboratory system. However, it seems reason- 
able to assume that the velocity of the a particle can 
be of the order of zero-point oscillation velocity of he- 
lium atoms, Av ~ hn 1 / 3 /m ~ 10 3 cm/s. The given 
estimate is based on the qualitative description [9( of 
the condensate wave function ^jy that varies substan- 
tially on the distances ~ n -1 / 3 . It means that the rel- 
ative uncertainty of the kinetic energy Em = Mv 2 a j2 of 
the a emitter with mass M should be not worse than 
AE M /E M ~ Av/v a ~ 10~ 6 . 

Obviously, we are dealing here with a factor that makes 
the detection of stimulated a emission difficult. The 
problem is not so much in the smallness of AEm / Em as 
in the nonconstancy of the kinetic energy of the moving 
nucleus in the superfluid helium (due to slowing down). 
One may hope, however, that this problem is solvable. 
Some ideas are proposed in the following, final section. 

CONCLUSION 

It is suggested that superfluid helium may influence 
the probability of nuclear a decay. Indeed, the super- 
fluid component consists of a significant part of helium 
atoms and their nuclei (a particles) in a single quantum 
state (in Bose condensate). Thus, the amplitude of a 
particle emission to this state is enhanced by the fac- 
tor of \JN C + 1, where N c is the macroscopic occupation 
number. 



It is assumed that the ion or neutral atom (with de- 
caying nucleus) moves through the layer of superfluid 
helium, so that the backward-emitted a particle is at 
rest (with respect to the condensate of 4 He atoms). It 
is shown that for a emitters with a lifetime 10~ 7 -10 -5 s 
the probability of stimulated a decay during the motion 
of the nucleus through the layer of thickness ~ 10 cm 
may be of the order of 10~ 3 -10 -2 . 

To detect the effect, it is necessary to compare prob- 
abilities of a decay for temperatures of helium above 
and below T\. The correlation would be an impressive 
demonstration of the possibility to influence a quantum 
process with characteristic energy release ~ 10 7 eV by a 
very weak impact at the level of T\ ~ 10~ 4 eV. 

Constant velocity of the moving a emitter (ion or 
atom) is necessary. Perhaps the energy losses of the ion 
due to slowing down can be recompensed by some ad- 
ditional ion acceleration directly in helium. Perhaps the 
ion or atom can be launched through the helium inside 
a nanotube without slowing down. Perhaps the ion or 
atom can be launched along a surface of superfluid he- 
lium, a small distance (of the order of 10 -8 cm) from 
the surface (in a bound quantum state corresponding to 
the movement normal to the surface). Thus, to detect 
the effect, significant efforts and additional ideas may be 
needed. 

It is useful to keep in mind that when studying 
new physical phenomenon something unexpected can be 
found. In particular, experimental study of the a decay 
of a nucleus at rest with respect to superfluid helium is 
of a certain interest (because the a particle leaving the 
Coulomb barrier is also "at rest" in some sense). 
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NS-3004.2008.2 and 2.1.1/4540 from MES of Russia. 
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